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Abstract 

The relativistic configuration interaction valence bond method has been used 

to calculate permanent and transition electric dipole moments of the KRb het- 

eronuclear molecule as a function of internuclear separation. The permanent 

dipole moment of the ground state X l Y> + potential is found to be 0.30(2) 

eao at the equilibrium internuclear separation with excess negative charge on 

the potassium atom. For the a 3 S + potential the dipole moment is an order 

of magnitude smaller (1 eao = 8.47835 10 -30 Cm) In addition, we calculate 

transition dipole moments between the two ground-state and excited-state 

potentials that dissociate to the K(4s)+Rb(5p) limits. Using this data we 

propose a way to produce singlet X l T, + KRb molecules by a two-photon 

Raman process starting from an ultracold mixture of doubly spin-polarized 

ground state K and Rb atoms. This Raman process is only allowed due to 

relativistic spin-orbit couplings and the absence of gerade/ungerade selection 

rules in heteronuclear dimers. 

PACS numbers: 03.75.Fi, 03.75.Be, 03.75.-b 
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In heteronuclear diatomic molecules the electrons can be distributed between the two 
nuclei unequally. As a result, there can be a small excess of negative charge near the more 
electronegative atom and thus heteronuclear dimers can have a permanent dipole moment. 
In this study we calculate the permanent dipole moments of the 1,3 E + electronic states of 
the KRb ground configuration, which dissociate to the K(4s)+Rb(5s) limit, and transition 
dipole moments between these states and excited states, which dissociate to K(4s)+Rb(5p) 
limits. Simple analysis of the ionization energies of K and Rb atoms, 4.339 eV and 4.176 
eV, respectively, shows that the K atom is most likely to have an excess of a negative charge 
when the two atoms interact. 

There are several applications of dipolar molecules in the physics of ultra-cold molecular 
gases. The long-range interaction between two dipolar molecules is governed by a electric 
dipole-dipole interaction and can, for example, significantly modify the many-body dynamics 
of trapped ultra-cold molecular Bose Einstein Condensates (BEC's) [1-6]. The atom-atom 
interactions in current atomic BEC's are spherically symmetric in origin. Following the 
proposal for homonuclear molecules of Ref. [7], dipolar molecules might also be formed in an 
optical lattice [8,9] An ultra-cold atom of each of the atomic species is held in an individual 
lattice site. A molecule in the electronic ground state can then be formed under the influence 
of laser light by inducing a two-photon transition. For any of the proposed applications, 
knowledge of the permanent or transition dipole moments is essential. 

The dipole moments have never been determined for KRb states of the ground configu- 
ration either theoretically or experimentally. We calculate the dipole moment of the singlet 
A Al £ + and triplet a 3 S + states of KRb as a function of internuclear separation R using an 
ab initio configuration interaction valence bond (VB) method [10]. Both relativistic and 
nonrelativistic calculations are presented, which allows us to determine the influence of rel- 
ativistic effects on the dipole moment. The use of configuration interaction (CI) in the VB 
method is essential due to the role of correlation in the formation of the molecular bond and 
nonzero dipole moment. The main contribution to the dipole moment of KRb is not due to 
the ground state atomic configuration, even though it predominantly determines the energy, 
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but is due to the population of excited valence and virtual orbitals. 

The accuracy of the calculation of the X X S + wave function is tested by comparing our 
ab initio potential with existing experimental [11-14] and theoretical [15-18] singlet poten- 
tials. We also present the a 3 £ + potential, compare it with the calculations of Refs. [17,18], 
and indicate the sensitivity of this shallow potential to the computational method. 

Furthermore, we analyze the possibility of the creation of dipolar KRb molecules in 
an optical lattice [8]. In this Paper we calculate the transition electric dipole moments 
between the X ' 3 S + states of the ground configuration and relativistic = ± , 1 components 
of excited 1,3 E + and 1,3 n states dissociating to the K(4s)+Rb(5p) limits, where Q is the 
projection of the total electronic angular momentum of the two atoms on the internuclear 
axis. This data provides us information about the most efficient scheme of forming KRb 
molecules in the singlet X X S + ground state from a K(4s)+Rb(5s) collision on the triplet 
a 3 E + potential. Moreover, we present our ab initio potentials for the excited states which can 
be used to identify the complex behavior of the transition dipole moment as function of R. 
Large-scale theoretical studies of the excited potentials of KRb were previously performed 
in Refs. [17,18]. It was shown in Refs. [19-22] that the dispersion interaction of excited 
KRb states is considerably stronger than in any other heteronuclear alkali dimer and that 
the KRb system is a good candidate for photoassociation experiments. Transition dipole 
moments to more-highly excited KRb states have been calculated in Ref. [17]. 

I. THEORY 

The electronic potentials and dipole moments are calculated with the configuration- 
interaction (CI) valence bond method. Both nonrelativistic and relativistic implementations 
of the method are used. The basic idea behind the valence-bond method is that the electronic 
molecular wavefunction |\I/ab) is constructed from wavef unctions, which are products of 
wavefunctions that describe the constituent atoms A and B. In essence, the molecular 
wavefunction is given by 
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\^AB)=Y. C «\ d < B ) (1) 

a 

where each \det AB ) is an antisymmetrized (A) product of two atomic Slater determinants 

\det AB ) = A(\det A ) ■ \det B )) . (2) 

The atomic Slater determinants \det A ) and \det B ) are constructed from one-electron func- 
tions centered on the nucleus of atom A and B, respectively. The variational CI coefficients 
C a in Eq. (1) are obtained by solving a generalized eigenvalue matrix problem, since in our 
approach the \det AB ) are not orthogonal. 

The permanent molecular dipole moment of the molecular wavefunction \^ab) is calcu- 
lated from 

N 

li= E Z k R k -e{^ AB \J2n\^A B ) (3) 

k=A,B i=l 

where N is the total number of electrons in the molecule and e is the electron charge. The 
first term of Eq. 3 depends on the position R k and charge Z k of nuclei k = A and B and 
is zero by choosing an appropriate coordinate origin. The second term of Eq. 3 depends 
on the electronic molecular wavefunction and the electron positions r*j. In the valence-bond 
method 

N 

{^ab\ £rl|^ B ) = NJ2 dffp AB (x,x) , (4) 

i=l K J 

where x = (r, k) denotes both the single-electron coordinate r and the variable k. For a 
nonrelativistic or relativistic calculation k has two or four values, respectively. The single- 
electron density matrix is 

p AB (x,x') = J2C:C^(x,x>), (5) 

a,/3 

where the transition density matrix between determinants \det AB ) and \det AB ) is defined by 

p a 'P(x,x') = (6) 

N 



D a p is the determinant of the N x N overlap matrix S"f = (ipi\ipj), and the ipi(x) = 
(x\ipi) and (pi(x) = (x\ipj) denote the single-electron functions used in constructing the basis 
functions \det^ B ) and \detp B ), respectively. Moreover, D aa and Dpp are the determinants 
of the matrices S"f = {4>i\4>j) and Sff = (<Pi\<Pj). The one-electron wavefunctions \ipi) 
and \ipj) are obtained by self-consistently solving Hartree-Fock and Sturmian equations 
for a nonrelativistic calculation and Dirac-Fock and Sturmian equations for a relativistic 
calculation [10]. Transition dipole moments can be derived in a similar fashion by assuming 
different initial and final states in Eq. 3. 



II. GROUND STATE DIPOLE MOMENTS 

The atomic determinants for the KRb dimer are constructed from single-electron Hartree- 
Fock or Dirac-Fock functions/orbitals for electrons in closed and valence shells. Sturmian 
wavefunctions complement the basis functions and are used to describe virtual orbitals. One- 
electron functions are characterized by a main quantum number n = 1,2, . . . and orbital 
angular momentum I = s,p, ... for both nonrelativistic and relativistic calculations. In 
addition for relativistic calculations the one-electron orbital is labeled by the total electron 
spin j. 

The closed shells Is 2 2s 2 2p 6 3s 2 of K and Is 2 2s 2 2p 6 3s 2 3p 6 3d 10 4s 2 of Rb form 
the core of the heteronuclear molecule. The superscript denotes the number of electrons in 
shell nl. For relativistic calculations 2p 6 is short for 2p\j 2 2p\j 2 etcetera. In our calculations 
excitations from these closed shells to valence and virtual orbitals are not included, i.e. all 
atomic determinants in the molecular basis contain the same number of electrons in these 
orbitals. Single electron excitations from the closed 3p 6 shell of K and 4p 6 shell of Rb 
are allowed and introduce core-valence correlations in the CI. The 4s, 4p, and 3d valence 
orbitals of K and 5s, bp, and id orbitals of Rb are allowed to contain at most two electrons. 
In addition, for the nonrelativistic calculation we allow at most two electrons in the Ad, bp, 
and 6p virtual orbitals of potassium and the bd, 6p, and 7p virtual orbitals of rubidium. For 
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the relativistic calculation computational limitations restricted us to the Ad and hp virtual 
orbitals of K and hd and 6p virtual orbitals of Rb. Both covalent and ionic configurations 
are constructed. The CI expansions with these basis sets have 493 nonrelativistic and 1459 
relativistic configurations. 

The ground state potentials have equilibrium distances R e and dissociation energies D e 
that agree with experimental RKR values [13,14] and other theoretical potentials. Figure 1 
shows the ground state X l Y, + and a 3 S + potentials of the KRb dimer as a function of the 
internuclear separation R. The solid curves of Fig. 1 describe results of our relativistic (Q 
= 0) ab initio calculation whereas the dashed curves are the theoretical data of Ref. [18]. 
The dotted line of Fig. 1 shows the singlet RKR potential of Ref. [14]. 

For the X X S + state the difference between the two theoretical calculations shown in 
Fig. 1 is about 3% at R e . A 20% difference exists at R e of the a 3 S + potential. Our estimate 
of the Cq/R 6 Van-der-Waals coefficient for the X : S + potential, obtained by fitting to the 
dispersion potential C 6 / R 6 + C 8 /R s + C 10 /R 10 ', agrees to 2-3% with the value of Derevianko 
et. al. [23]. The experimentally determined dissociation energy for the X 1 E + potential 
[13,14] is 26 cm -1 smaller than ours. 

We performed calculations of the X X E + and a 3 £ + state potentials and their permanent 
dipole moments in both nonrelativistic and relativistic approximation. For the same set of 
one-electron orbitals the influence of relativistic effects on the potential energy is very small 
whereas its influence on the dipole moments is large. At R e a relativistic calculation leads 
to a 30% increase in the absolute value of the dipole moment for the X l Y> + state and a 
50% decrease for the a 3 £ + state. Moreover, we find large correlation effects in the electric 
dipole moment in both nonrelativistic and relativistic calculations. The implementation of 
core polarization, for instance, leads a 50% of reduction in the absolute value of the dipole 
moment at the bottom of the singlet potential. 

Figures 2 and 3 show the electric dipole moments of the X 1 S + and a 3 S + potentials 
of KRb as a function of internuclear separation, respectively. A negative dipole moment 
implies an excess electron charge on the potassium atom. The presented dipole moments 
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are calculated with the relativistic and nonrelativistic Hamiltonian. Curve 1 in Figs. 2 and 
3 is determined from the most accurate nonrelativistic calculation, which includes the two- 
electron occupation of all orbitals upto 6p for the K atom and 7p or the Rb atom. Curve 
2 in Figs. 2 and 3 shows the f2 = relativistic calculation. The number of excited orbitals 
is smaller than in the nonrelativistic basis set and limited to excitations upto 4d and hp for 
K atom and 5d and 6p for Rb atom. The difference between the dipole moments of the 
Q = and 1 component of the a 3 S + potential is much smaller than the uncertainties of our 
calculations. Tables I and II tabulate the permanent dipole moment of the X X E + and a 3 S + 
potentials shown in the figures. 

Our calculation shows that the distribution of the charge density between K and Rb 
is very diffuse. It means that the dipole moment depends not only on the electron charge 
transfer from Rb to K atom but also on the induced polarization from the charge transfer, i.e. 
the electrons occupy excited p and d orbitals. Charge transfer and the induced polarization 
are of opposite sign. This is especially true for the a 3 £ state. We find that the excess charge 
is almost equal for the singlet and triplet states, but the latter has a smaller dipole moment. 

The permanent dipole moments of the X X S + and a 3 £ + state multiplied by R 7 , based on 
a the relativistic calculation are shown in Fig. 4. It shows that the long-range behavior of 
the dipole moment of the two states are equal and proportional to 1/R 7 . This long-range 
behavior is due to the modification of the molecular wave function by the dipole-dipole and 
dipole-quadrupole multipole interactions [24]. 

We believe that our convergence with respect to correlation or the number of basis func- 
tion is of the same order as the difference between nonrelativistic and relativistic calculation. 
Consequently, we feel that the most accurate dipole moment is obtained from averaging the 
nonrelativistic and relativistic calculation and assuming an one-standard deviation uncer- 
tainty equal to half the difference between the two calculations. For the X 1 E + state the 
dipole moment is —0.30(2) ea at R e = 7.7 a , while for the a 3 E + state the dipole moment 
is -0.02(1) ea at R e =ll.2 a . 
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III. TRANSITION DIPOLE MOMENTS 



In this section we analyze the possibility of the creation of dipolar KRb molecules in an 
optical lattice. We assume that colliding atoms are initially in the doubly spin-polarized 
state, which only allow them to come together on the state a 3 S + potential. Thus photoasso- 
ciation via an excited molecular 3 S or 3 I1 state produces a 3 S + molecular levels. In principle 
this makes the production of molecules in the X X S + state problematical, since one has to 
get from the triplet to singlet spin manifold. There is, however, a viable route from the dou- 
bly spin-polarized colliding ground state atoms to ground X X S levels via the excited state. 
When the detuning of the photoassociation laser from atomic resonance is small compared 
to the excited state spin-orbit splitting, as will be the case for most practical photoassocia- 
tion schemes, then an excited O = or 1 Hund's case (c) state can have both singlet and 
triplet character. Thus, an excited molecular state can be formed from excitation from a 
a 3 E + state that can re-emit light in a transition to a X l lu + state. This process is absent 
in homonuclear dimers, since the additional gerade- ungerade selection rule for electronic 
transitions prevents it. 

Figure 5 shows the Q = O 1 * 1 and 1 relativistic excited state potentials as a function of 
internuclear separation. The excited state potentials are obtained with the basis described in 
Sec. II. At short internuclear separation the potentials are described by the ( 2S+1 )A ± Hund's 
case (a) labeling following Ref. [17] where A is a projection of the electronic orbital angular 
momentum on the molecular axis and S is the total electron spin. At longer R the curves are 
better described by fi ± Hund's case (c) labeling and relativistic effects are important. The 
potentials dissociate to the excited K(4s)+Rb(5p 2 Pj) and K(4p 2 Pj )+Rb(5s) fine structure 
limits. The two lowest dissociation limits correspond to the two fine structure states of the 
excited Rb atom plus a ground state K atom. The long-range behavior of the relativistic 
potentials dissociating to these two limits is attractive. There are seven attractive potentials 
with f2 = and 1. Our potentials agree with the calculations of Refs. [17,18]. 

The attractive excited-state potentials are the most likely candidates for use in two-color 
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Raman photoassociation experiments. We calculate transition dipole moments relevant for 
transitions from the a 3 £ + (0~, 1) state, through the attractive f2 = ± and 1 potentials, to 
the X 1 E + (0 + ) state. Figures 6 and 7 show the nonzero transition dipole moments from the 
a 3 S + and X : S + states, respectively. Curves with the same style in the two figures indicate 
the same intermediate excited state. Photon selection rules ensure that + — > CT transition 
are not allowed. At long-range the transition dipole moments become independent of R and 
their absolute values approach the Rb 5s — > 5p( 2 Pj) transition dipole moment when the 
corresponding excited potential dissociates to the K(4s)+Rb(5p 2 Pj) limit. 

At short-range internuclear separation the dipole moments strongly depend on R. This 
behavior reflects the change from a Hund's case (a) to a relativistic coupling scheme between 
20 a and 30 a . In Fig. 6 the a 3 S+(fi = 1) to 2 1 S+(0+) and in Fig. 7 the X 1 S+(0+) 
to 2 3 S + (1) transition dipole moments approach zero at short R because singlet to triplet 
transitions are not allowed. In Fig. 6 both a 3 S+(l) to 2 3 S+(1) and a 3 S+(0") to 2 3 S+(1) 
have a small dipole moment at small internuclear separation. This can be understood by 
noting that the structure of the ground and excited state potentials shown in Figs. 1 and 5 
is similar to that of homonuclear alkali-metal dimers, where gerade/ungerade symmetry is 
valid. In homonuclear dimers electric dipole transitions between gerade and ungerade states 
are forbidden. 

The sudden change in dipole moment in Fig. 7 near 13 a is related to the avoided 
crossing indicated by the circle in Fig. 5 between the Q = 1 components of 1 3 II and 2 3 E + 
potentials. We estimate that the uncertainty of the transition dipole moments is 0.1 ea 
based on a comparison of the calculated dipole moments at R = 100 ao with the known 5s 
to 5p dipole moments of Rb [25]. 

IV. CONCLUSION 

We determined the permanent dipole moments of the 1,3 S + states of the ground configu- 
ration of the KRb heteronuclear molecule using a nonrelativistic and relativistic CI valence 



9 



bond method. The KRb permanent dipole moments are small compared to "truly" polar 
molecules, such as NaCl. It might, however, be large enough for experiments that aim to 
confine KRb in optical traps. In addition we calculated the potential energy curves and 
transition dipole moments to excited states correlating to K(4s)+Rb(5p) atomic limits. We 
have shown that there exist allowed transitions starting from colliding doubly polarized K 
and Rb atoms via an exited state to the singlet X 1 E + state. 

The calculation of the electric dipole moments is a first step towards obtaining quanti- 
tative estimates of photoabsorption and molecular production rates in a gas of K and Rb 
atoms. In the future we plan to evaluate Frank-Condon factors between vibrational lev- 
els of ground and excited state potentials. In addition we need to investigate the effect of 
black-body radiation on the dipolar molecule by evaluating Frank-Condon factors between 
vibrational levels of the ground state. 
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TABLES 

TABLE I. The electric dipole moment of the state. Both nonrelativistic (nrel) and 

relativistic (rel) results are presented. The internuclear separation R is in units of ao and the 
dipole moment is in units of eaQ. 



R 


nrel 


rel 


6 


-2.52x 10 _1 


-2.82x 10 1 


6.5 


-2.60x 10 _1 


-2.89x 10 _1 


7 


-2.69x 10 _1 


-3.00x 10 _1 


7.5 


-2.77x 10 _1 


-3.12x 10 _1 


8 


-2.81x 10 _1 


-3.22x 10 _1 


8.5 


-2.80x 10 _1 


-3.28x 10 _1 


9 


-2.72x 10 _1 


-3.28x 10 _1 


9.5 


-2.57x 10 _1 


-3.20x 10 _1 


10 


-2.35x 10 _1 


-3.02x 10 1 


11 


-1.77x 10 _1 


-2.45x 10 _1 


12 


-1.18x 10 _1 


-1.75x 10 _1 


13 


-7.25x 10~ 2 


-1.13x 10 _1 


14 


-4.25x 10~ 2 


-6.80x 10 2 


15 


-2.61x 10~ 2 


-4.11 x 10~ 2 


16 


-1.56x 10~ 2 


-2.36x 10~ 2 


18 


-6.44 x 10~ 3 


-9.44 x 10- 3 


20 


-2.48x 10~ 3 


-2.48x 10 -3 


22 


-LlOx 10~ 3 


-LlOx 10- 3 


24 


-5.49x 10" 4 


-5.49x 10~ 4 


26 


-3.01 x 10" 4 


-3.01 x 10" 4 


28 


-1.74x 10~ 4 


-1.74x 10~ 4 


30 


-1.08x 10~ 4 


-1.08x 10~ 4 
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32 -6.88x 1(T 5 -6.88x 1(T 5 

34 -4.57x 1(T 5 -4.57x 1(T 5 

36 -3.02x 1(T 5 -3.02x 1(T 5 

38 -2.06x 1(T 5 -2.06x 1(T 5 

40 -1.43x 1(T 5 -1.43x 1(T 5 
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TABLE II. Electric dipole moment of the a 3 E + state. Both nonrelativistic (nrel) and rela- 
tivists (rel) Q =0 results are presented. 



R nrel rel 

5 -2.22x 1(T 2 -1.20x 1CT 1 
5.5 -3.94x 10~ 3 -1.13x 

6 1.66x 10~ 2 -1.05x 10" 1 
6.5 3.21x 10" 2 -8.21x 10~ 2 

7 3.94x 10~ 2 -5.66x 10~ 2 
7.5 3.90x 10~ 2 -4.54x 10~ 2 

8 3.32x 10~ 2 -4.19x 10~ 2 
8.5 2.46x 10~ 2 -4.05x 10~ 2 

9 1.56x 10~ 2 -3.91x 10~ 2 
9.5 7.26x 10~ 3 -3.71x 10~ 2 

10 4.50x 10~ 4 -3.45x 10~ 2 
10.5 -4.67x 10~ 3 -3.16x 10~ 2 

11 -8.18x 10 -3 -2.86x 10~ 2 
11.5 -1.03x 10~ 2 -2.56x 10~ 2 

12 -1.14x 10~ 2 -2.28x 10~ 2 
12.5 -1.17x 10~ 2 -2.01x 10~ 2 

13 -1.13x 10~ 2 -1.77x 10~ 2 
13.5 -1.07x 10~ 2 -1.54x 10~ 2 

14 -9.80x 10~ 3 -1.34x 10~ 2 
14.5 -8.82x 10~ 3 -1.16x 10~ 2 

15 -7.83x 10~ 3 -9.02x 10~ 3 
18 -3.24x 10~ 3 -3.24x 10~ 3 
21 -1.22x 10~ 3 -1.22x 10~ 3 
24 -4.78x 10~ 4 -4.78x 10~ 4 
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27 -2.04x 1(T 4 -2.04x 1(T 4 

30 -9.67x 1(T 5 -9.67x 1(T 5 

33 -5.18x 1(T 5 -5.18x 1(T 5 

36 -2.72x 1(T 5 -2.72x 1(T 5 

39 -1.54x 1(T 5 -1.54x 1(T 5 

42 -9.14x 1(T 6 -9.14x 1(T 6 
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R (units of a Q ) 

FIG. 1. The ground state and a 3 S + potentials of KRb as a function of internuclear 

separation. Solid curves are obtained in the present study, dashed lines are the results of Ref. [18], 
and the dotted line is the RKR potential of Ref. [14]. (1 ao=0.0529 nm) 
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R (units of a Q ) 



FIG. 2. The electric dipole moment of the X l Y, + ground state of the KRb dimer as a function 
of internuclear separation. The dipole moment is found from a nonrelativistic (curve 1) and a 
relativistic (curve 2) calculation. 
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5 10 15 20 

R (units of a Q ) 

FIG. 3. The electric dipole moment of the a 3 S + O = state of KRb as a function of internuclear 
separation. The dipole moment is found from a nonrelativistic (curve 1) and relativistic (curve 2) 
calculation. 
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FIG. 4. The relativistic electric dipole moment multiplied by R 7 of the X 1 T, + and a 3 £ + (f2 = 0) 
states of KRb as a function of internuclear separation. 
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FIG. 5. Electronic SI = ± , 1 potentials of excited states of the KRb dimer in an intermediate 
region of internuclear separation. The short-range n^ 2S+1 ^A^ and long-range m(n ± ) labeling is 
indicated. The numbers n and m show the energy-ordered appearance of these states. 
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FIG. 6. Transition dipole moments between the ground a 3 E + and excited states of the KRb 
dimer as a function of internuclear separation. The curves are labeled by short- and long-range 
symmetries as in Fig. 5 
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FIG. 7. Transition dipole moments between the ground X 1 E + and excited states of the KRb 
dimer as a function of internuclear separation. The curves are labeled by short- and long-range 
symmetries as in Fig. 5 
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